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ABSTRACT 

High energy emission from blazars is thought to arise in a relativistic jet launched by a supermassive 
black hole. The emission site must be far from the hole and the jet relativistic, in order to avoid 
absorption of the photons. In extreme cases, rapid variability of the emission suggests that structures 
of length-scale smaller than the gravitational radius of the central black hole are imprinted on the jet as 
it is launched, and modulate the radiation released after it has been accelerated to high Lorentz factor. 
We propose a mechanism which can account for the acceleration of the jet, and for the rapid variability 
of the radiation, based on the propagation characteristics of large-amplitude waves in charge-starved, 
polar jets. Using a two- fluid (e^) description, we find the outflows exhibit a delayed acceleration 
phase, that starts when the inertia associated with the wave currents becomes important. The fluids 
propagate with the wave at approximately the sonic speed, corresponding to a bulk Lorentz factor 
7 « \{)'^/S.t\^f^Kr^^'^ L^J^ Mg^^^ out to radius ri « /S.t\^f^K^^^ l]^^ M^^^ ^c, after which the Lorentz 
factor accelerates as 7 oc r. (A^ioo is the variability time in units of 100 s, the pair multiplicity 
at one gravitational radius, the "47r-luminosity" of the jet in units of 10^^ erg/s, and Mg the 
black-hole mass in units of 10^ Mq.) The time-structure imprinted on the jet at launch modulates 

photons produced by the accelerating jet provided K^g < 14 Atiooilg^Mg"^, suggesting that very rapid 
variability is confined to sources in which the electromagnetic cascade in the black-hole magnetosphere 
is not prolific. 

Subject headings: MHD - plasmas - waves - BL Lacertae objects: individual (PKS 2155-304) - galax- 
ies: jets - gamma-rays: galaxies 
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1. INTRODUCTION 

Observations by the H.E.S.S. collaboration of TeV 
gamma-ray emission from the blazar PKS 2155-304 
reveal very rapid variability fAharonian et al.l 120071 : 
IHESS Collaboration et alll2010( ). at very high fiux lev- 
els. In the most extreme flare, variations on a timescale 
of a few hundred seconds at a fiux level correspond- 
ing to an isotropic luminosity of 10"*^ erg/s were mea- 
sured. If, as expected, the mass M of the central 
black hole is 2 x 10^ M©, these observations imply 
structure in the jet that is roughly one hundred times 
smaller than the gra vitational radius = GM/c^ 
(|Begelman et al.ll2008f) . Though this is the most extreme 
example, several other blazars exhibit ver y rapid vari- 
ability in GeV and TeV gamm a-rays fe.g.. lAlbert et al.) 
I2007t lAckermann et al.l l20ldl ). which it is proving dif- 
ficult to accommodate in the standard synchrotron- 
self-Compton picture, mainly because of the very high 
Lore ntz factor and l ow magnetisation required of the 
iet (iLevinsonl 120071: iBoutelier et all 120081: iGraff et al " 



2008 ; IKatarzvnski et al.l 120081: IKusunose &: Takahara 



20081: 'Mastichiadis fc Moraitis"2008': 'Neronov et al .ll200 
GhiscUini et al. 2009; Giannios et al. 2009; Pag gi et al.l 
20091: iTammi fc Duffy 120091: iRieger fc VoW 120101: 
Nalewaiko et al.. 2010, ') 



In the framework of ideal MHD, radial (uncollimated) 
relativistic jets do not accelerate after they pass through 
the fast magnetosonic point. Collima tion, however, re- 
quires special boundary conditions (jLvubarskvl I2009L 

|john.kirk@mpi-hd.mpg.de, iwona.mochol@mpi-hd.mpg.de| 



120101 ). It is, therefore, difficult to envisage the production 
of a jet with high Lorentz factor and low magnetisation. 
An isolated, impulsive, ejection event in a n ideal MHD 
flow is able to circumvent this problem (jGranot et al.l 
I20T0I) . but, in a fluctuating jet, dispersion filters out 
the small timesca le structure, and limits the acceleration 
(|Levinsonll2010al ). It appears, therefore, that it may be 
necessary to go beyond the ideal MHD approximation in 
order to understand the observations of very rapid vari- 
ability in blazars. 

Non-ideal MHD effects can become important when 
plasma is charge- starved: If the number density of 
charged particles is limited, this places a maximum on 
the absolute value of the charge density, which can, 
for example, result in the inability of the plasma to 
screen out the component of the electric field parallel 
to the magnetic field. Charge starvation also places 
an upper limit on the available current density. As 
this limit is approached, the relative drift-speed of the 
charged components becomes relativistic, and their as- 
sociated inertia begins to contribute to the fiuid stress- 
energy tensor. This situation might plausibly arise 
in a black-hole magnetosphere, because axisymmetric, 
general-relativistic, MHD simulat ions of jets launched 
by accreting, rotating black holes (|De Villiers fc Hawlevl 
120031: iMcKinnev fc Gammij I2004D reveal a conical re- 
gion around the rotation axis into which the accreting 
plasma does not penetrate. As in a pulsar magneto- 
sphere, the matter density in this region is likely to be 
determined not by accretion, but by the rate at which 
electron-positron pairs are created in the strong elec- 



tromagnetic fields that penetrat e it (iGoldreich fc Julianl 
[1961 IBlandford fc ZnaiekllTOTTl: ILevinsod I2000ir " If an 

outflow results, the density of pairs decreases, and, far 
from the hole, non-MHD effects connected with particle 
inertia can become important. 

In the case of an axisymmetric, force-free magne- 
tosphere, it is known that plasma is ejected, carry- 
ing off energy mainly in the form of Poynting flux via 
the "Blandford-Znajek" mechanism. On the axis itself, 
the energy flux vanishes, so that the polar regions of 
the jet do not dominate the overall energetics. How- 
ever, observations of rapid variability suggest that ax- 
isymmetry may not be a good approximation, since 
they imply small-scale structure in the black-hole mag- 
netosphere. Simulations in which the black-hole spin 
and the asymptotic magnetic held are misaligned find 
a Poynting fl ux compa rable to the Blandford-Znajek 
value (Palenz uela et al.i 2010), making it plausible that 
also more complex non-axisymmetric held structures can 
power a substantial, magnetically dominated, polar jet. 

In the following, we develop this idea by examining the 
propagation characteristics of nonlinear electromagnetic 
waves above the polar regions of a rotating black hole. 
Using a model consisting of cold electron and positron 
fluids, we show that a circularly polarised magnetic shear 
propagates radially outwards at roughly the fast magne- 
tosonic speed until it reaches the point where the ideal 
MHD description loses its validity. In (non-ideal) MHD 
language, this happens when the inertia of the plasma 
particles begins to affect the conductivity, i.e., when the 
drift sp eeds implied by the plasn ia current become rela- 
tivistic (|Melatos fc Mehoselfl996HMeieiil2004[ ). This can 
occur at a large distance from the black hole, depend- 
ing on the density of injected pairs. The wave then goes 
through a phase of delayed acceleration in which it con- 
verts Poynting flux into kinetic energy flux. If particles 
radiate gamma-rays in the acceleration zone, then, de- 
spite the large spatial extent of the source, the spatio- 
temporal structure of the shear wave, that is imprinted 
on it close to the black hole, modulates the radiation, 
provided the mass-loading of the jet is sufficiently small. 

In section [2] we discuss the parameters used to spec- 
ify the physical conditions in the jet. The two-fluid jet 
model is presented in section [31 First, the nonlinear 
plane- wave solution representing a magnetic shear is de- 
rived, then radial propagation in spherical geometry is 
discussed. It is shown in the Appendix that general rela- 
tivistic effects drop out in the short-wavelength approx- 
imation (c/oj <C r) when the Kerr metric is used. The 
application to blazar variability is discussed in sectional 
and our conclusions summarised in [5] 

2. JET PARAMETERS 

Consider a radial outflow consisting of an electron- 
positron plasma emerging from the polar regions of the 
magnetosphere of a rotating black hole, and denote by 
L and M the total luminosity and mass-flux carried in 
a solid opening angle fls- The physical conditions in the 
flow can be specifled via three dimensionless parameters: 

1. The nonlinearity or strength parameter a is a di- 
mensionless measure of the energy-flux density. 
For a circularly polarized vacuum electromagnetic 
wave of frequency u!/2tt and electric field E, the 



strength parameter is conventionally defined as 
a = eE/{mcoj), and measures the Lorentz factor 
that an electron would achieve if it were acceler- 
ated from rest over a distance of (27r) times one 
wavelength in the field E. Assuming radial propa- 
gation, the corresponding luminosity is 



(1) 



and we use this expression to define a for a general 
(non-vacuum) wave. In the absence of radiation 
losses, a cx r"^, and is determined by specifying 
its value at some fiducial radius. For this we chose 
ro = c/w, although our treatment is, of course, 
valid only for r 3> ?'o and we certainly do not expect 
radial flow to extend to such small radii. With this 
choice, flQ is independent of uj: 

(2) 



a = ao (ro/r) 
Go -■ 



-,1/2 



3.4 X 10"L^/' , 



(3) 



where L46 = (47r/f7s) L/ (lO"^ erg/s) is the 
"isotropic" or "47r" luminosity of the jet, scaled 
appropriately. The energy radiated per unit solid 
angle by a jet is directly measurable if the distance 
to the object is known. For the rapidly variable 
gamma-ray flare of PKS 2155-304, it corresponds 
to an isotropic luminosity of roughly lO'*^ erg/s, so 
that, for this object, L46 > 1. 

2. The mass- loading of the wind is conventionally de- 
scribed by the /n-parameter introduced by iMicheil 
(fl969l) : 



fi = L/Mc^ 



(4) 



In the case of an electron-positron jet, /i denotes 
the Lorentz factor each particle would have if the 
entire luminosity was carried by a cold, unmagne- 
tised flow. It is constant in those parts of the jet 
in which pair creation and radiation losses can be 
neglected. 

The magnetisation parameter a describes the ra- 
tio of the energy flux carried by electromagnetic 
fields to that carried by particles. For monoener- 
getic electrons and positrons of Lorentz factor 7, 



Wt) - 1 



(5) 



In a cold, non-accelerating, ideal MHD flow, a 
is constant (assuming pair creation and radiation 
losses are negligible). However, as we show below, 
a is not constant in charge- starved jets, even in the 
absence of dissipation. For this reason, we specify 
its value at the "launching" radius, inside of which 
the ideal MHD approximation is assumed to hold, 
and denote this quantity by ctq, even though the 
region of constant a is unlikely to extend to radii 
as small as r ~ c/w. The particle Lorentz factor at 
the launching point is then 



70 = m/ (c^o + 1) 



(6) 
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The mass-loading parameter fj, is determined by the 
physics of the pair-production cascade close to the black 
hole. A more intuitive measure, therefore, is the pair 
multiplicity k, which relates the pair (proper) num- 
ber density n± to the number density of electrons (or 
positrons) needed to screen out the (magnetic-)field- 
aligned component of the rotation-induced electric field. 
Adopting the definition conventionally used in pulsar 
physics, but replacing the a ngular velocity of the n eu- 
tron star by c/vg gives (e.g., iLvubarskv feKiilllffll 



'-l±n± 



B 



2'Ker„ 



(7) 



where ^± is the Lorentz factor of the fluids. In the inner 
regions of the flows we consider, where ct = cto S> 1, 
the fluids move non-relativistically in the wave frame, so 
that K « a/ (4/i). Thus, in the absence of radiation losses 
and pair production, k cx r^^ in this region. Physically, 
it is the value of k at the outer boundary of the pair- 
production region that is most relevant. This is thought 
to be close to the black hole, but its precise location is 
unknown. In the following, therefore, we specify k by its 
value Kr at r = r„: 



c 

Lor. 



(8) 



3. THE TWO-FLUID MODEL 



The simplest model of an electron-positron plasma 
that captures the physics connected with the finite in- 
ertia of the charge-carriers is that of two cold, oppo- 
sitely charged fluids (denoted by suffices — and -|-). 
We ad opt this niodel , embed it in a Kerr metric, fol- 
lowing [Khanni (119981) and Koide (2009), and look for 
large-amplitude waves propagating in the radial direc- 
tion in Boyer-Lindquist coordinates. To keep the analy- 
sis tractable, only transverse waves with vanishing phase- 
averaged components of the electric and magnetic fields 
are treated. In this case, since V • £^ = 0, the number 
density is the same for each fluid. We further assume 
the wave carries no radial current, so that the radial 
fluid velocities also equal each other: Vf+ = Vf- = Vf, 
and restrict the treatment to waves in which the merid- 
ional and azimuthal fluid velocities are equal in magni- 
tude but of opposite sign: Vg^ — —Vg_, v^_^ — —v^_. 
It follows that the fluids have the same Lorentz factor: 

/ \ 1/2 

7+ 7- = 7 = c/ f - Vf - vj^ - w?^ j , and the 

same radial component of the four-velocity, which we 
write as a dimensionless momentum: p|| = Vf"f/c. Cir- 
cularly polarised waves are likely to be the most impor- 
tant in the polar regions of a rotating black hole, and 
these are best treated by introducing complex quantities 
to describe the transverse components of the fluid mo- 
menta: pi_ = (vg^ + iv^^'j 7/c (= - (vg^ + 7/c) 
and the electric and magnetic fields: E = Eg + iE^, 
B = Bg + iB^. In Appendix El we derive the continu- 
ity equation, the equations of motion of the fiuids and 
the two relevant Maxwell equations (Faraday's law and 
Ampere's law) in the small- wavelength approx imation, 
r ^ c /u, starting from the formulation given bv lKhannal 
(|1998D . In the lowest order, we search for large-amplitude 



plane- wave solution s using the approach introduced by 
lAkhiezer fc PolovinI (|1956[ ). 



3.1. Nonlinear waves 

Expressing all quantities in terms of the phase (f> de- 
fined in (jASP, t he continuity equation (|A12[) and Fara- 
day's law (|A13p integrate immediately to give 



nA = constant 
B = iE/l3^ , 



(9) 
(10) 



where A = 7 — (p||//3w), and c/3w is the phase speed of 
the wave. The equations of motion to this order, (jAlsp 
and (|A16p . are 



ujA'^ = — —lm{p^B* 
d0 mc 

dpj_ . eB 
ujA—— = —i — A . 

dffl mc 



(11) 
(12) 



Solutions to these equations can be found with super- 
luminal phase speed, /3w > 1, bu t these wa ves do not 
propagate close to the black hole (|Kirk|[2010[ ). Here, we 
concentrate on subluminal waves, for which the condition 
A = holds. Physically, this means that the particles are 
in resonance with the wave, i.e., the radial components of 
the fluid velocities equal the phase velocity of the wave. 
In this case (|9]) shows that the phase dependence of the 
density is arbitrary, (IT^ is trivially satisfled, and (|lip 
requires that the transverse components of the fluid ve- 
locities are parallel to the magnetic field: Im {p±B*) = 0. 
Thus, the plasma current is directed along the magnetic 
field and, to this order, the forces exerted on the fiuids by 
the fields vanish. Finally, the current and magnetic field 
are linked by Ampere's equation, which, to this order, 
reads 



dB 8T:inecp±/3^ 



^(1-/3^) 



(13) 



Combined with (llip this implies that the magnitude of 
the magnetic field is phase-independent. Viewed from a 
frame that moves radially with speed /3w , the electric field 
vanishes, and the wave is simply a static magnetic field 
of constant magnitude whose direction rotates through 
27r radians over one wavelength. At each point, the cur- 
rent and, hence p±, is parallel to the magnetic field to 
zeroth order. The rate at which the B-vector rotates is 
arbitrary, being determined by the dependence of the 
fluid density n on phase. In the following, we select 
the simplest case, where n, \B\^ and \p±\^ are all con- 
stant and the wave is a monochromatic magnetic shear: 
B (xp±(x e±**. 

3.2. Radial evolution of a magnetic shear 

The slow evolution of the subluminal magnetic shear 
wave as it propagates outwards at r 3> Tg is governed by 
the flrst-order equations in the expansion in e ~ c/ (ujr), 
as derived in Appendix |Al Making obvious simplifica- 
tions to the notation, these are the continuity equation 
(IM) : 



P|| /R 



(14) 



4 




4 5 6 7 
log(R) 

Fig. 1. — The magnetisation parameter cr, transverse fluid mo- 
mentum px and Lorentz factor of the wave, 7w, as functions of the 
dimensionless radius R = firuj/ (aoc), for ap = 10^ and fi = 10^^ . 
The approximate solutions given in Eq. II20I I are also shown. 



the equation of energy flux conservation: 

^ = 7(1 + cr) 

and the radial momentum equation 
di^ R\p±\^ 



(15) 



dR 



(16) 



where the momentum-flux density per unit mass-flux 



v=P\\ 1 



and Lj = Lj/uJp (with the plasma frequency defined us- 



ing the proper fluid density: 



Sirne /m), R 



(lor/c) (/i/ao) is the radius in units of the critical radius, 
inside of wh ich the superluminal modes do not propagate 
(|Kiii[20T0l ). Note that in a non-monochromatic wave, 
the quantities p\\, w, cr and 7 are replaced in these equa- 
tions by their phase-averages. 

According to ([T^ . tr, as defined in (|A32p is related to 
the fluid momentum components through 



(18) 



-1/2 



where 7w = (l — The condition that the wave 

velocity equals the radial component of the fluid velocity, 
A = 0, implies 



7w = 7' 



7(i + Kr 



(19) 



The five equations p^ - ([T^ . together with the def- 

inition 7 = (l+pn +|p^| j , determine the radial 

dependence of the six unknown wave variables 7w, 7, p\\, 
\p±\, w, and cr. It is straightforward to reduce these to a 
first-order ordinary differential equation for p±(7w), for 
example. Solutions extend from i? = to i? = 00 pro- 
vided they are launched at super-magnetosonic speed: 
= (a^/t) - 1 < M^^^- At i? -J> 0, \p±_\^Q and 7 7w. 
whereas at i? —>■ cx), cr and 7 — ?> /x, so that the wave 
converts all of the Poynting fiux to kinetic energy at large 
radius. 

The radial dependence of the wave goes through three 
phases. At small R, \p±\ <^ 1 and the wave is essentially 



a cold MHD structure in which the inertia associated 
with the current is negligible. There is no acceleration of 
either the wave speed or the fluids in this regime, and the 
magnetisation parameter a remains constant at its initial 
value ctq. Assuming ctq 3> 1, this region is restricted to 
R <^ fi/ao. At intermediate radii, one readily finds an 
approximate solution: 



p±^l 
7w2i? 

valid in the range 



(20) 



(21) 

Finally, at large radius, i? ^ /i, only kinetic energy re- 
mains: 7 ~ 7w ~ \p±\ ^ 1- This behaviour is illus- 
trated in Fig. [TJ 

4. APPLICATION TO BLAZAR VARIABILITY 

The locations of the three phases of wave propagation 
illustrated in Fig. [T] depend on the parameters oq, fJ, (or 
f^rg): Co and uj. As discussed in Sect. [H it is possible 
to infer values for oq and w directly from the observed 
flux and variability timescale. Another parameter may 
be eliminated by fixing the wave speed at its launching 
point. A very slow, sub-magnetosonic outflow can be 
described by the for ce-free MHD e quations, and would 
accelerate as 7 oc r (|Bucklevlll977[ ). until it approached 



(17) the sonic speed, where 70 « cr. 



1/2 



On the other 



hand, all waves launched at super-magnetosonic speeds 
(erg < /i^/"^) behave similarly, as described in Sect. 13. 2[ 
with the acceleration phase moving out to larger radius 
as the initial magnetisation decreases. It suffices, there- 
fore, to analyse the case of mildly supersonic launch: 
corresponding to the maximum magnetisa- 



2/3 



tion of a supersonic flow. 

However, the uncertainty associated with the unknown 
mass-loading of the flow can be removed only by mod- 
elling the pair cascade. Leaving this quantity as a free 
parameter, the radius race at which acceleration begins, 
corresponding to i? sa m/'^'Oi is 



1/3 2/3 / I \-l/3 

-1 2Af^/\2/3rl/6,.2/3 
' 1.^ ^EjqqK 1^4g ilia 



pc 



(22) 



where we define the variation timescale in units of 100 s 
to be Aioo = (27r/tL') / (100 s), and write the mass of the 
black hole as M = Mg x 10^ M©. For r < race the Lorentz 
factor of the fiow remains constant at roughly the sonic 
speed: 



70^6.5 X lO^AtK^ln-^/^ 



1/3 



(23) 



and, for r > Ta 
with r: 



the Lorentz factors increase linearly 



7w«7/V2 « 7.4 X 10^ {r/lpc)K;^Mg^ 



(24) 



The solutions presented in section 13.21 eventually con- 
vert all of the Poynting flux to kinetic energy flux at 
large radius. However, in the case of blazars, this is un- 
likely to be realised, since the resulting Lorentz factor 
(= fj,) is very large. Instead, dissipative processes so far 
neglected, such as instabilities in the wave-solution, or 
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interaction of the jet with the external medium, or with 
ambient photons, are hkely to intervene. 

The wave propagates radiaUy with fixed frequency. 
However, if it converts part of its energy into high- 
frequency (3> Lu), forwardly beamed photons via an emis- 
sivity that is modulated by the wave-phase, then the dif- 
ference between photon and wave propagation speeds will 
lead to a smoothing of the modulation in the observed 
photon signal. This loss of short-timescale variability 
becomes more effective as the size of the radiating sec- 
tion of the wave increases. Similarly, if the photons do 
not propagate exactly in the radial direction, smoothing 
of the modulation will be produced by the difference in 
light-travel time to the observer from different parts of 
the spherical wavefront. It is straightforward to derive a 
criterion on the size of the emitting region (assumed ~ r) 
and the Lorentz factor of the jet, such that fluctuations 
of frequency y are not suppressed in the p hoton signal 
(|Miche]|ll97lllAronsllT979l: IKirk et al.ll200l : 

-fl2TTc/Lu>r . (25) 

In the acceleration region, 7^ oc r, so that this condi- 
tion is fulfilled everywhere within this region, provided it 
is satisfied at the beginning, where r — race, and 7w = 7o- 
Combining (l22t and p3|) the requirement that modula- 
tion on a timescalc of lOOAiioo seconds should not be 
filtered out leads to an upper limit on the multiplicity: 

K,^<14Aiiooi^f Afg-i . (26) 

Equation (j26p implies that electron-positron pair cre- 
ation is much less effective in the central engine of a 
rapidly vari able blazar than it is in a pulsar magneto- 
sphere (e.g.. lMedin fc Laill2010( ). but this is perhaps not 
unexpected, given that a neutron star surface is able to 
anchor a very strong magnetic field. However, it also 
implies that blazars exhibiting extreme variability con- 
tain a charge -starved magnetosphere able to support a 
vacuum gap (iLevinsonI 120001 l2010bl ). This scenario is 
particularly attractive, because the non-stationary na- 
ture of gap di scharges found in puls ar-related studies 
(|Levinson et al. 2005; TimokhinI I2010D , suggests a nat- 
ural source of short-timescale (< rg/c) variability in the 
outflow from a black-hole magnetosphere. 

5. SUMMARY AND CONCLUSIONS 



In this paper, we describe a mechanism that causes a 
magnetically dominated, radial outflow from a black-hole 
magnetosphere to enter a delayed acceleration phase, 
starting at a distance from the hole given by (|22|) . Ap- 
plying this mechanism to blazar jets, we derive a con- 
straint, (|26|) . on the pair density in the magnetosphere 
that would allow radiation produced where the jet ac- 
celerates to retain any short-timescale structure imposed 
on it close to the launching site. 

The mechanism is based on an analysis of the prop- 
agation characteristics of a nonlinear wave - specifi- 
cally a circularly polarised magnetic shear ~ in a low- 
density plasma. Such a wave, we suggest, is likely to 
be launched in the polar regions of a rotating, accret- 
ing black hole, and, in a non-axisymmetric picture, may 
fluctuate on a time shorter than r^/c, as indicated by ob- 
servations of the source PKS 2155-304. Acceleration is a 
result of charge-starvation - a non-MHD effect that arises 
when the relative drift-speed of the oppositely charged 
constituents in a low-density plasma becomes relativis- 
tic. The analysis employs a cold two-fluid model of the 
plasma, and uses a short- wavelength perturbation expan- 
sion to find the evolution of the radially propagating, 
nonlinear wave. The equations are derived in Kerr geom- 
etry. However, under the conditions we envisage, where 
the wavelength of the oscillation is of the same order 
in the expansion parameter as the gravitational radius, 
general relativistic effects do not appear in the governing 
equations. 

Several important problems remain to be investigated. 
These include the nature of the dissipation and radia- 
tion mechanisms, and the effect these might have on the 
propagation of the wave, as well as the possibility of 
modelling the multi-wavelength blazar spectrum. Fur- 
thermore, although the picture of a circularly polarised 
magnetic shear that is static in the jet frame is intuitively 
attractive, this is only one specific, nonlinear solution of 
the governing equations; other polarisations and other 
modes, such as the line arly polarised "striped wind" 
(jLv ubarskv fc Kirl3 120011 ) or the electro magnetic mode 
of superluminal phase-speed (|Kirkl'2010f) may also prove 
important. Nevertheless, the underlying physical cause 
of the acceleration — the inertia of the charge-carriers — 
suggests that delayed jet-acceleration may be a generic 
phenomenon. 



APPENDIX 

EQUATIONS OF WAVE PROPAGATION IN KERR GEOMETRY 

Here we derive the equations governing the radial propagation of transverse, circularly polarised, electromagnetic 
waves in a plasma consisting of cold electron and positron fluids that are embedded in a Kerr metric. We use a 
short-wavelength approximation: c/wr ~ e <C 1, where uj is the wave frequency, and assume the gravitational radius 
is of the same order in e as the wavelength. In this appendix we s et G = c = 1, s o that with M the black-hole mass, 
M ~ l/o; ~ er. We start from the two-fluid equations as given by iKhannal ()1998[ ). and use the (essentially standard) 
notation of that paper. 

As measured by a fiducial observer corotating with a black hole, the fluid 3- velocities and Lorentz factors, and the 
electric and magnetic fields are denoted by v±, 7±, and E and B respectively, the suffices -|- and — denoting the 
positron and electron components (charge ±e, mass to). The continuity equation (|Khannalll998l Eq. (22)) reads 

{dt-13- V) n±7± + V • {an±-i±v±) = , (Al) 

where V is an operator in curved 3-dim space described by Qij, a is the lapse function and /3 the gravitomagnetic 



6 

potential. The equations of momentum conservation for each fluid (jKhannalll998l . Eq. (23)) 

a-^ {dt- l3-V)S±=eg + ^ ■ S±- a'^V ■ (aV"±) 

±e"/±n±E ± ej±n±v± x B (A2) 

reduce to the equations of motion 



a 



dt 



±- (E' + e'"'vj±Bk) (A3) 



in the special case of cold, collisionless fluids (in Khanna's notation e — j±n±m, S± = jj_n±mv± and ^ = 
n±'mjj_v±v±). Faraday's law and Ampere's law are 

V x{aE)=-{dt-C(3)B (A4) 
y X (aB) = {dt - Cp) E + inae {n+p+ - n-p-) , (A5) 

where Cf3 is a Lie derivative along (3: 

CpE ^{(3-V)E-{E-V)(3 . (A6) 

We now restrict the treatment to radially propagating waves by assuming the field and fluid variables depend only 
on t and r, and introduce the wave phase a function of t and r that depends on the (as yet unspecified) wave phase 
velocity /3^{r), which is a function of r alone. The phase of an outwardly propagating vacuum wave in Kerr geometry 
is ^vac = w (i — r*), where 

rr „/2 I „2 

:dr' , (A7) 



2Mr' + a- 



a (in this appendix) is the Kerr parameter, and uj is the wave frequency measured by an observer at infinity (e.g. 
iThorne fc Pricdll986l Eq (8.66). In analogy with this expression we write the phase of the nonlinear wave as 



Because we restrict the treatment to transverse waves, E and B are automatically divergence-free, and the plasma 
is charge-neutral. The waves of interest have Vf+ = Vf- = Vf and Vg^ — —Vg_, v^^ — ~v^_, so that 7+ = 7- = 7, 
n+ = n_ = ri, and we describe them using complex quantities for the transverse (dimensionless) momenta p± — 
7 (vg^ + iv^^^ and for the electric and magnetic fields: E = Eg + iE^, B = Bg + iB^. In accordance with standard 
notation, we rename the radial component p|| — pf. 

Transforming the independent variables in (jAlHASp from (t, r) into a "fast" phase variable and a slow radial coor- 
dinate: (0, p), where p = er, we now expand in the small parameter e, assuming p ^ r^. Keeping terms of zeroth and 
first order, the derivatives are replaced according to 

d d 
dt-^^d-^ 
d dp d d<j) d 
dr dr dp dr d4> 

d / 2M\ UJ d 

-%-['^'-)^;d-^ ^^''^ 

Expanding the dependent variables according to p\\ = p''^^ +^?'||^^ '^^^ finds the zeroth-order equations are those 
of continuity: 

d 



^£(„(o)A(o))=0 , (A12) 
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Faraday's and Ampere's laws: 



ui 



/3w 



(A13) 



/3w ^'Z' 

and momentum/energy conservation: 

o (0) 

^AWZL + Aim fpf =0 (A15) 

O0 m \ / 

o (0) 

^^(0)^ „ A f7(")i?(°) =0 (A16) 

dcp m \ II / 

^A(°)^ -^Refpf =0 , (A17) 

where A — 7 — p||//3w The monochromatic, subluminal sohition to these equations has A^"-* — 0, and \B''^^\ , , 

(0) " 
P± 



2 

, pfi"' and n'"^ all independent of ^ 



I 

Taking account of this, the first-order equation of continuity is: 



Faraday's and Ampere's laws are: 



90 p i9p 



/3w 90 



i|.(pB(0)) „ .^SvrenWpf =0 , (A20) 
pop P 



and momentum/energy equations give: 



o (1) 



d(j) m \ 



''dp m p V / p 



a (1) Q (0) 

d<h 00 m \ " " ) 



,(0)^ + A^ f^(o)^(o) +,J0)^(0)A - zcote^^ (A22) 

" op m p \ II / p p 



A(0)^ + A(i)^ - ^Re fp(°)£;W%pW£;(")*) 
9© 00 m \ J 



d(j) dcj) m 

+Pr^ + --Re(pri?(°)V0 . (A23) 
II 9p m p \ J 

The "slow" dependence of the zcroth-order quantities on p follows by elimin ating secular terms in the first-order 
quantities, i.e. by imposing the condition that they are periodic in 0. Equation (|A18P can immediately be integrated 
over 0, yielding, when periodicity is imposed, 

lArp2^(oUo)Uo , (A24) 



p2 dp " ^1 

Similarly, (jA19[) integrates to give 
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However, this merely constrains the average components of the wave fields, which we assume to vanish. In order 
to integrate (|A23p and (IA2ip . it is first necessary to use Ampere's law (|A20I) to re-express p^j^'' in the expressions 
Im (p'f^B'-"^*^ and Re (pf^E'^^l*^ in terms of 5*°) and respectively. One then finds 



d_ 

dp 



2 / (0) (0) (0) , 

\ II 8nm 



= 



d_ 
d'p 



IQiTtn 



(A26) 



(A27) 



Equations (|A24p . (jA26|) and (jA27p suffice to determine the dependence on p of the phase-averaged, zeroth-order 
variables. Note that, to this order, M does not appear; i.e., general relativistic effect s do not enter. Furthermore, 
because we assume cold, dissipationless fiuids that interact only via the wave fields, (|A26I) simply states the conservation 
of the sum of the zeroth-order particle and field contributions to the phase-averaged energy fiux in fiat space, expressed 
in differential form. 
Integrating (|X24| and (|M6)) . 



and from (|A27p 



where we define 



2mr^n(")p[|°^=M/a 



(0) 



1 



2Pl 



(A28) 
(A29) 

(A30) 
(A31) 



M/flg is the mass-flux per unit solid-angle, p — L/M \s the mass-loading parameter, v is the radial momentum flux 
density per unit rest mass and the magnetisation parameter, defined as the ratio of the field and particle terms in the 
energy fiux density, is 



|S(0)/^(0)| 
(0) 

777, 



(A32) 



Using the definition of the strength parameter ([3]) enables the mass-fiux to be expressed in terms of ap and /i, leading 
to 



P 



(0). 



(A33) 



1/2 

where Wp = (Snne'^ /rnj is the "proper" plasma frequency. 
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